Theory of the field- induced BEC in the frustrated spin-| dimer compound BaCuSi 2 6 
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Building on recent neutron and NMR experiments, we investigate the field-induced exotic criticality observed 
in the frustrated spin-| dimer compound BaCuSi2 0e using a frustrated model with two types of bilayers. A 
semiclassical treatment of the effective hard-core boson model shows that perfect interlayer frustration leads 
to a 2D-like critical exponent <f> — 1 without logarithmic corrections and to a 3D low temperature phase with 
different but non vanishing triplet populations in both types of bilayers. These further suggest a simple phe- 
nomenology in terms of a field-dependent transverse coupling in the context of which we reproduce the entire 
field-temperature phase diagram with Quantum Monte Carlo simulations. 
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Introduction — Field-induced Bose-Einstein condensation 
(BEC) of the excitations in antiferromagnetic (AF) spin sys- 
tems is a fascinating phenomenon. Theoretically predicted 
some years ago [1, 2], it was first observed experimentally 
in the spin-i dimer compound TICUCI3 [3, 4, 5]. Recently, 
a very intense activity has emerged [6], showing a large and 
rich variety of field-induced bosonic phases in such quantum 
materials. Interestingly, a quantitative theoretical description 
is often possible through effective theories written in term 
of hard-core bosons. Such theories have been quite success- 
ful to explain various non-trivial features such as magnetiza- 
tion plateaus in spin ladders [7, 8] and frustrated lattices like 
Sr2Cu(BC<3)2 [9, 10], field-induced Luttinger liquid behavior 
in quasi-lD systems [11], or even a more exotic phase, analog 
to a supersolid, where a uniform triplet condensate coexists 
with a spontaneous spatial modulation of the triplets [12]. 

Another remarkable behavior has been reported in 
BaCuSi206, a 3D stack of spin-^ bilayers with frustrated 
interlayer coupling known as Han purple. This compound 
undergoes a field-induced AF ordering described as a triplet 
BEC [13], but the field dependence of the critical temperature 
above the quantum critical point (QCP) at H c ~ 23.2 T has 
an anomalous exponent: [14] 

Tbec(^) ~ (H - H c f (1) 

with cf> ~ 1. Since the BEC has an exponent = 2/D, 
this behavior is typical of a 2D system. This lead Sebas- 
tian et al. [14] to interpret their result as a frustration driven 
"dimensional reduction" of the interlayer coupling jj_ based 
on the observation that the triplet transverse dispersion ~ 
jj_ cos (k x /2) cos (ky/2) cos k z is flat at the AF condensate 
wave vector (71", 7r), so that there is no 3D coherence of the 
condensate. This was later confirmed by inelastic neutron 
scattering experiments [15], which indeed reported an al- 
most fiat 3D dispersion at (tt, tt), but which also revealed 
a richer structure with inequivalent bilayers of at least two 
types: J A = 4.27(1) meV, J B = 4.72(1) meV [15], as de- 
picted in Fig. 1(a). The presence of at least two types of lay- 
ers has been confirmed by high field NMR experiments which 
have identified a larger triplet population in A-layers than in B 
ones. [17]. These NMR results point however to the 3D nature 



of the ordered phase at low temperature since B-layers have a 
small but non zero population above H c . 

Taken together, these experimental results [14, 15, 17] lead 
to a fundamental puzzle: How can a 3D bosonic population 
lead to a 2D critical exponent for the BEC transition line? 
And more generally, what is the exact role played by the frus- 
tration? Two theories have been put forward, with somehow 
incompatible conclusions. Batista et al. [16] have investigated 
a model with equivalent layers. As ancitipated in Ref. [14], 
they found that, although the low temperature phase is 3D, the 
critical temperature does not have the 3D exponent <p = 2/3. 
They obtained <p — 1, but with logarithmic corrections. This 
conclusion holds provided the interlayer coupling is fully frus- 
trated. By contrast, Rosch and Vojta [18] have investigated 
a model with different bilayers and came to the conclusion 
that the interlayer coupling cannot be fully frustrated since 
this would be inconsistent with the NMR observation of a non 
zero triplet population in B-layers just above H c . 

In this Letter, we show that all these experimental data can 
in fact be reconciled in the context of a model including the 
two types of bilayers and perfect frustration. In particular, 
we show that the presence of two types of bilayers lead to 
Tbec(H) ~ (H—H c ) without logarithmic corrections, while 
the B-layers get a triplet population immediately above H c . 
However, this population only grows as pB oc (H — H c ) 2 , i.e. 
much more slowly than that in the A-layers (pa oc H — H c ). 
Starting form a 3D frustrated spin bilayer model with two 
types of bilayers, these conclusions have been reached by in- 
vestigating the spin-wave (SW) corrections around the classi- 
cal ground-state (GS) of an effective frustrated bosonic model 
for the triplets with an energy barrier between the A and B 
planes of A ~ 0.45 meV [15]. These results can be rein- 
terpreted as a field-dependent effective quantum tunnelling in 
the transverse direction ~ (H — H c ) for the triplets. This 
naturally implies a 3D critical BEC temperature with an ef- 
fective 2D-like exponent <f> = -g + | and has allowed us to 
reproduce the entire experimental T — H phase boundary us- 
ing large scale quantum Monte Carlo (QMC) simulations of a 
3D effective bosonic model with a field-dependent tunneling. 
Effective model — The theoretical analysis starts from a mag- 
netic Hamiltonian describing a 3D array of frustrated bilayers 
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with two types of dimers [Fig. 1(a)] that we can decompose 
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with J A ~ 4.27 meV, J B 
meV forBaCuSi 2 6 [15]. 
governed by 
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The frustrated inter-bilayer part is 
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Note that the space coordinates are labelled by r, I where I = 
1 , 2 are the layer indices within each bilayers, and elementary 
displacements r and S are defined in Fig. 1(f). Since we con- 
sider the strong coupling limit (Ja,b ^ ■/")> eacn bilayer can 
be described by a low energy theory where only the two lowest 
states of a dimer (the singlet \s) = (|tl) — lit)) A/2 an d the 
polarized triplet = Iff) shown in Fig. 1(d)) are retained, 
each dimer being replaced by a single site with hard-core bo- 
son degrees of freedom (|0) = |s) and |1) = \t^}) [2, 7]. Each 
AF bilayer is then described by interacting hard-core bosons 
on a 2D square lattice, and the 3D effective frustrated bosonic 
Hamiltonian [Fig. 1(b)] reads H bos =Ha + Hb+H±, with 
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where /i^ = J a — 77 and hb = Jb — H = ha + A.. The 
frustrating inter-layer part mixes bosons a and £>: 
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Mean-Field analysis — We first analyse the problem using a 
mean-field (MF) Hartree-Fock approach, searching for a GS 
wave function of the form "to = i{ja <8> V'b with |?/a) = 
IIjga ( sin ^1°) i + eiVj cos 8 -f I ■ When decoupled, each 
layer displays successive phases upon changing the chemical 
potential. If /i < —2J" (> 2jH) the bosonic system is empty 
(full): = 7T (0 = 0) and |$ = ^ |0)j (\i/>) - Uj 
whereas for — 2j" < /x < 2j" the system is compress- 
ible and the bosons can form a T = Bose-Einstein con- 
densate with < 9 < ir and a fixed (pinned) phase ip (fi- 
nite superfluid stiffness). Interestingly, the layers A and B 
start to accommodate bosons for different values of [i, mean- 
ing different values of the magnetic field 77, resulting in the 




FIG. 1: Schematic models for BaCuSi20e: (a) Frustrated bilayer 
array with two types of dimers; (b) and (c) Effective bosonic models 
in the field-induced critical regime: (b) Frustrated hard-core bosons 
model with different chemical potential; (c) Non-frustrated hard-core 
boson model with a field-dependent effective transverse hopping (see 
text); (d) The 4 energy states of isolated dimers A and B versus the 
external magnetic field; (e) Triplet dispersions in isolated bilayers 
Ea = 7" (cosk x + COSA;,,) and 75b = Ea + A induced by the 
inter-dimer couplings j"; (f)View of the frustrated lattice (a) from 
above with the unit vectors f and 8. Bottom: Angle representation 
of the classical GS. (g) fj, < fi c : both A and B layers are empty, (h) 
fi c < ft < /i c + A: A layers start to accomodate bosons while B 
layers are still gapped. 



distinct triplet branches shown in Fig. 1(e). When frustra- 
tion [Eq. (6)] is taken into account, the gap remains open at 
the MF level because frustration leads to a cancellation of 
the hopping processes between A and B and only allows for 
non-zero diagonal interactions which just renormalize the gap 

A = ( J B - J A ) / U - • Therefore, at the MF level the 
triplet densities pa(b) = (1 + cos #a(b))/2 appear at different 
critical fields with a quite sizable gap A ~ 3.5 T. However, 
as already discussed by Rosch and Vojta [18], the existence 
of a second critical field where triplets would start to populate 
the B planes is ruled out by NMR results [17] which clearly 
show a unique critical field with a finite, while reduced, triplet 
population on the B layers. Therefore, one needs to go beyond 
MF and consider the effect of quantum fluctuations around the 
classical GS in order to understand how the triplets populate 
the B sites, and how they propagate along the 3D direction. 
Linear spin-waves — We replace a hard-core bo- 
son by a spin-1/2 with classical coordinates S — 
1/2 (cosy sin 9, sin <p sin 9, cos 9), leading to the vector 
representation depicted in Fig. l(g-h). We develop the SW 
calculation around these classical vectors, including the finite 
external field [19, 20]. After a few steps, we compute the 
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SW-corrected densities p A w and p| w which both appear 
for p > fi c . While p A w remains linear with p — p c , SW 
fluctuations close the "classical gap" A and induce a small 
triplet population in B layers as soon as p > p c . This bosonic 
occupation being classically forbidden in the region (b) 
p < fi c + A, it is purely induced by quantum fluctuations. 
More precisely, an unexpected quadratic dependence on the 
external field is found, with the following dependence 



p| w ~ const x 



A 



(a*-Mc) 2 



(7) 



k = (6J") _1 being the compressibility of an isolated layer, 
and const is 0(1) (see Fig. 2). SW results are shown in 
Fig. 2(a) with realistic parameters for BaCuSi2C>6 (J\\ = 0.5 
meV and A = 0.45 meV) and various values of the frus- 
trating coupling Jj_. Such a quadratic dependence, which is 
compatible with current NMR data [17, 21], has a dramatic 
effect on the 3D criticality of the system. Indeed, consider- 
ing the non-frustrated version of this problem: a system hav- 
ing a similar energy barrier A between adjacent planes but 
with non-frustrating inter-layer tunnelling tj_, one obtains a 
unique QCP with a B occupation growing linearly like /?b — 
t]_K,(iJ, — fi c )/ A 2 . This implies that the quadratic population 
in the frustrated case Eq. (7) results from an effective field- 
dependent inter-layer tunnelling amplitude t c f ~ y/p — p c . 
Therefore the 3D coherence, governed by the 2 nd order hop- 
ping process [Fig. 1(c)], is controlled by the energy scale 
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const x (J±) K 2 (/i — p c ). 
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As an alternative check for such a field-dependence for 
the effective 3D tunnelling amplitude, we introduce a Peierls 
phase twist ip± in the transverse direction and compute the 
transverse superfluid stiffness = d 2 (Tt(p±)) / dip 2 ^ | _ Q . 
A finite stiffness signals the appearance of 3D phase co- 
herence, hence of an effective hopping i 3 o with ~ 
i3Dsin 2 9a- SW results for — /3^/[pa(1 — Pa)] are shown 
in Fig. 2(b) for realistic parameters. The agreement with the 
linear behavior Eq. (8) is excellent. On the other hand, equiv- 
alent layers (A = 0) do not lead to a simple linear 3D hopping 
but to quite strong logarithmic corrections (inset of Fig. 2). 
Consequences for the BEC critical temperature — The kinetic 
energy of a diluted isotropic Bose gas is e(q) = f^q 2 / (2m*). 
If D> 2 there is a finite-T BEC when the De Broglie length 
becomes greater than the inter-particle distance d ~ P^, 
at T c oc ■ In our specific case, close to the critical 

field, the kinetic energy is anisotropic: e(q) = t (q 2 + q 2 ) + 
t3D? 2 . The linear t3o(p) leads to an effective mass m* ~ 
(p — p c )~n, thus implying a modified 3D BEC temperature 
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with an effective exponent <f> = + — 1, as experimentally 
measured [14]. It is important to note that with only one type 
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FIG. 2: Scaling plots from SW calculations with realistic parameters 
Jy = 0.5 meV, A = 0.45 meV and various values of the frustrat- 
ing couplings indicated on the graph, (a) Quadratic behavior of 
p| W Eq- (7); (b) Linear behavior for the 3D tunnelling t3D Eq. (8) 
deduced from the transverse stiffness px- Black lines are fits using 
Eqs. (7,8) with const = 1.6. Inset: Comparison between equivalent 
layers (A = 0, closed symbols) and non-equivalent layers (A = 0.5 
meV, open symbols). The absence of finite size effects and the linear 
behavior t^o oc p — fi c in the latter case contrast with the strong 
logarithmic corrections of the former case. 



of layer (A = 0), the logarithmic corrections discussed above 
for t3D will show up in T C (H), as also found in Refs. [16] 
using a different theoretical approach. 

Quantum Monte Carlo — A direct test of these conclusions 
on the model of Eqs. (2,3) with QMC is impossible because 
of the severe minus sign problem induced by frustration. In 
order to test the internal logic of the previous arguments re- 
garding the linear field-dependence of the critical temperature 
Eq. (9), we performed QMC simulations on the non-frustrated 
effective bosonic model [Fig. 1(c)] HeS = "Ha + Ti-B + 7~t±, 
where Ha and Hb are given by Eqs. (4) and (5), and Ti ± is 
the non-frustrated version of Eq. (6). Taking realistic param- 
eters for BaCuSi 2 6 : J 11 = 0.5 meV and A = 0.45 meV, 
the transverse hopping is adjusted such that the entire ex- 
perimental phase boundary is reproduced. The inverse map- 
ping bosons — > dimers is done such that the T = critical 
field H c takes the experimental value ~ 23. 2T for = 0. 
We achieved large-scale QMC calculations on systems of size 



L x L x L/2 (8 < L < 48) with < tf < 0.05 meV. 
The BEC transition was revealed by stiffness crossings of 
p s (T c ) x L, as expected for a 3D-XY transition. Shown 
in Fig. 3(a-b), the experimental data points are reproduced 
for the entire BEC critical dome using a continuously vary- 
ing transverse coupling t c f plotted in Fig 3(c). A fit to 
Eq. (8) leads to a rough estimate of the frustrating coupling 
Jj_ ~ 0.15(4) meV. As a comparison, the 2D BKT transi- 
tion for a single layer with J" — 0.5 meV is also displayed 
in Fig. 3(a). The curvature is clearly incompatible with the 
experimental data, thus demonstrating that the linear behavior 
of T c cannot be simply attributed to a 2D BKT transition. 
Summary and experimental consequences — We have shown 
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FIG. 3: (a) Temperature-magnetic field phase boundary for 
BaCuSi20e. Experimental data [14] (open circles) are confronted 
to QMC results (colored squares) obtained with a field-dependent 
transverse hopping t±. The solid line is a guide to the eyes for the 
linear part. The 2D-BKT transition line of a single decoupled bi- 
layer is also shown for comparison, (b) Zoom over the critical linear 
regime where more QMC data points are shown, (c) Effective hop- 
ping t c f used to fit the experimental phase boundary. The red curve 
is the theoretical t c ± using Eq. (8), the best fit being obtained for 
He ^ 23.5(3) and jf ~ 0.15(4) meV. 



that the combined effect of two types of bilayers and of inter- 
layer frustration in BaCuSi2C>6 solves all puzzles regarding 
the physics close to the critical field H c . It explains the linear 
behavior of T c with H — H c reported in Ref. [14], removing 
the logarithmic corrections of the single-bilayer model inves- 
tigated in Ref. [16], and it accounts for a small but non-zero 
triplet population in the B-layers. These results also point to 
a phenomenological description in terms of a field-dependent 
effective transverse hopping < 3 rj ~ H — H c which reconciles 
the concept of "dimensional reduction" at the QCP (<3D = 
at H c ) with the 3D coherence of the low-T phase above H c . 

This model also makes a number of predictions of experi- 
mental relevance. Firstly, the triplet population in the B-layers 
is quadratic with H — H c . This behavior is compatible with 
current NMR data [21], but more precise NMR measurements 
are needed to actually check this power law. Secondly, our es- 
timate w 0.15 meV implies for the transverse dispersion at 
(k x ,k y ) = (0,0) a band-width of 2(j£) 2 /A « 0.1 me V, that 
can be measured with inelastic neutron scattering. Finally, this 
model allows in principle to calculate the T = properties at 
high field, i.e. beyond the field where classically the B-layers 
get populated. A proper treatment of the vicinity of this field 
requires a priori to go beyond linear SW theory and is left 
for future investigation. It would be very useful however to 
measure the magnetization at very low temperature to see if a 
trace of this classical transition can be detected. 
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